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ABSTRACT 

We study the radio-far infrared (FIR) correlation in “blue cloud” galaxies chosen from the PRism 
MUltiobject Survey (PRIMUS) up to redshift (z) of 1.2 in the XMM-LSS field. We use rest-frame 
emission at 1.4 GHz in the radio and both monochromatic (at 70/im) and bolometric (between 8 — 

1000 /im) emission in the FIR. To probe the nature of the correlation up to z ~ 1.2, where direct 
detection of blue star-forming galaxies is impossible with current technology, we employ the technique 
of image stacking at 0.325 and 1.4 GHz in the radio and in six infrared bands, viz. 24, 70, 160, 250, 

350 and 500 /rm. For comparison, we also study the correlation for more luminous galaxies that are 
directly detected. The stacking analysis allows us to probe the radio-FIR correlation for galaxies 
that are up to 2 orders of magnitude fainter than the ones detected directly. The fc—correction 
in the infrared wavebands is obtained by fitting the observed spectral energy distribution (SED) 
with a composite mid-IR power law and a single temperature greybody model. We find that the 
radio luminosity at 1.4 GHz (Li.4ghz) is strongly correlated with monochromatic FIR luminosity 
at 70 /im (Lyo^jm) having slope 1.09 ± 0.05 and with bolometric luminosity (Ltir) having slope 
1.11 ± 0.04. The quantity (/tir(= logio[LTiR/(3.75 x 10 ^^Li.4ghz)]) is observed to decrease with 
redshift as (/tir oc (1 -I- probably caused due to the non-linear slope of the radio-FIR 

correlation. Within the uncertainties of our measurement and the limitations of our flux-limited and 
color-selected sample, we do not find any evolution of the radio-FIR correlation with redshift. 

Subject headings: methods : data analysis - techniques : image processing - surveys - dust - galaxies 
: ISM - galaxies : statistics - infrared : galaxies - radio continuum : galaxies 


1. INTRODUCTION 

The radio-far infrared (FIR) correlation is one of the 
tightest observed correlations in astrophysics that con¬ 
nects several independent physical parameters in the in¬ 
terstellar medium (ISM). The radio luminosity and the 
FIR luminosity of star-forming galaxies are observed to 
be cor related over five orders of magnitude for the global 
scale (iHelou et abl 119851: lCondonl[i992t lYun et al.l 120011 : 

lAppleton et al.l 12004 ISargent et al.l I2010H with disoer- 
sion less than a factor of 2. The radio luminosity is typi¬ 
cally measured at 1.4 GHz and the FIR luminosity can be 
both monochromatic (at 24, 60 or 70/nn) or bolometric 
(between 40 and 120/im or between 8 and 1000/nn). 

The radio-FIR correlation is well studied for galaxies 
in the local universe for several classes of galaxy mor¬ 
phology like spirals, ellipticals, dwarf irregulys, etc. It 
is known to hold good at global dCondonlllOf^lYun et al.l 

l200H[Price fc Duridll992HWunderlich et al.lll987ll as well 

as at local scales (few 100 pc to few kpc) within galax¬ 
ies (iBasu et al.l 1201 2bt iDuinas et aH 120111 : iHughes et al.l 
120061 : iHoernes et al.lll998[l . At the brightest end of FIR 
luminosity, the relationship is observed to hold for (ul¬ 
tra) luminous infrared galaxies [(U)LIRG] and star-burst 
galaxies. At the faintest end it holds in dwarf galax ies 
(|ChvW et al.ll20lit iRovchowdhurv fc Chengaluill2012ll . 

It is believed that star-formation connects the two 
regimes of emission. Synchrotron (also referred to as 
non-thermal) emission in the radio band is caused by 


acceleration of cosmic ray electrons (CREs) in the galac¬ 
tic magnetic field produced by supernova explosions of 
massive stars. In the FIR, the emission originates due 
to dust re-radiation, heated by ultraviolet (UV) photons 
from massive (> 10 Mq), short lived (^ 10® yrs) stars. 
However, the tightness seen in the correlation needs to 
be explained, as a number of independent physical quan¬ 
tities are responsible for the emission in each regime 
like, the magnetic field, number density of GREs, energy 
losses of GREs, star formation history, dust/gas den¬ 
sity, dust absorption efficiency, etc. Several models have 
been proposed to explai n the tightness seen in the radio- 
FIR correlation (see e.g.: rV51klll989HHelou fc Bicavlll993t 
iNiklas fc Beckl[ 997fl. More recent models bv iLacki et al l 
(|2010t) and iBelll (|2003f ) have shown that the above men¬ 
tioned factors conspire to maintain the tightness ob¬ 
served for the global radio-FIR correlation. 

Observationally, it is important to assess the form of 
the radio-FIR correlation at high redshifts as it might de¬ 
pend on the evolution of ISM parameters with redshift 
{z) like synchrotron and inverse-Compton losses, dust 
content, star formation rate, magnetic field strength and 
overall SEP (see e.g.: iMur phvi l2009t iLacki et ahl l201d: 

iSchleicher fc Becfl l20iab Recently, iSchleicher fc BecS 

(1201311 predicted a modification of the form of the radio- 
FIR correlation, based on the observed relationship be¬ 
tween magnetic field strength and star formation rate 
caused due to turbulent amplification of the magnetic 
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field. A breakdown in the correlation is expected depend¬ 
ing on the dominant energy loss mechanism of the CREs 
in the radio domain, i.e., synchrotron, inverse-Compton, 
bremsstrahlung and/or ionization losses. 

Typical (Icr) sensitivity of most of the existing deep ra- 
dio surveys are lim ited only to fe w tens of uJv (see e.g., 
iBondi et al.l 120031 [VLA-VVDS ] : iSchinnerer et al.1120101 
[VLA-COSMOS]; iHodge et al.l 120111 fEVLA-StriDe82]. 
etc.). However, a few deeper surve ys exists reaching 
Icr sensitiv ity < 10 /rJy (seee.g., iMiller et al.l 120131 
[E-CDFS]; iMorrison et al.l 120101 [GOODS-N]). These 
obseryations can detect normal galaxies (Ti, 4 GHz 
10^2 w Hz“^) up to redshift of ~ 0.2 at 1.4 GHz 
with > 5(7 sensitiyity, making it difficult to study the 
radio-FlR correlation for such galaxies at higher red- 
shifts. The correlation has been studied for (U)LlRGs 
with higher lumi nosity iL^ 4 GHz 10^^ W Hz~^) up t(3 
redshifts of ^ 3 (lApoleton et al.l 12004 iMao et ahl 120111 : 
iDel Moro et al.ll2013ll . Such galaxies can haye significant 
contamination due to AGNs and compact nuclear star- 
bursts. Eyen in the case of relatiyely low optical depth, 
starburst related free -free absorption can giye rise to sub¬ 
stantial obscuration (jGondon et ahlFToMIl that can affect 
the form of the correlation. It is therefore imperatiye 
to study the radio-FlR correlation for less extreme star¬ 
forming galaxies at high redshifts where the bulk of the 
radio and FIR emission originates from star formation. 

In this paper, we study the properties of the radio-FIR 
correlation, both the slope and the traditionally defined 
‘g’ parameter, for a flux limited and color selected sam¬ 
ple in the XMM-LSS field. We explore the correlation 
for blue star-forming galaxies up to z ~ 1.2 employing 
the technique of image stacking. Due to the inherent 
flux limitation of the parent sample, we detect normal 
star-forming galaxies up to z ~ 0.9 and more luminous 
galaxies aboye that. For comparison, we study the cor¬ 
relation for luminous galaxies that are directly detected 
in this field up to z ^ 0.95. 

The paper is organized as follows: In Section 2, we 
describe our sample selection and data. We discuss the 
technique of image stacking at 0.325 GHz and 1.4 GHz 
in the radio and at 24, 70, 160, 250, 350 and 500 fxm 
in the FIR and the fc—correction method in Section 3. 
We present our results in Section 4 and discuss them 
in Section 5. Throughout this paper, we assume a flat 
A—GDM model with Hq = 70 km s“^, Dm = 0.27 and 
Da = 0.73. 

2. SAMPLE SELECTION AND DATA 

To study the radio-FIR correlation at high redshift, 
one requires robust identification of galaxies along with 
accurate spectroscopic redshifts. The radio emission 
in galaxies mainly originates from non-thermal emis¬ 
sion and thermal bremsstrahlung emission. The thermal 
emission can affect the form of the radio- FIR correlation 
(|Hoernes et al.l[T998t iHughes et al.ll2006ll . Howeyer, ow¬ 
ing to the steep spectrum, more than 95 percent of the 
total emission at low r adio frequencies 0.325 GHz) is 
non-thermal in origin (jBasu et al.ll2012ah . Thus, to en¬ 
sure that the bulk of the radio emission is non-thermal in 
origin, a deep radio survey at low frequencies, like 0.325 
GHz, is necessary. Finally, a deep survey in the FIR 
regime at longer wavelengths (A > 20 fim) is important 
to avoid contamination from polyaromatic hydrocarbon 


features in the 5 — 10 fim regime up to z = 1. 

Here, we combine our Giant Metrewave Radio Tele¬ 
scope (GMRT) observations at 0.325 GHz with archival 
mid-infrared, far-infrared and 1.4 GHz radio data in the 
XMM-LSS field (centered at RA=2h 21m 00s, Dec=—4° 
30' 00" J2000). In Table [D we present the salient fea¬ 
tures of the various multiwavelength surveys used in this 
paper and Figure [T] shows the footprints of these surveys 
overlaid on the 0.325 GHz GMRT mosaic image. 

2.1. PRIMUS galaxy sample 

2.1.1. Salient features of the PRIMUS survey 

Our parent galaxy sample is d rawn from th e 
PRIsm MUlti-object Survey iPRIMUSi lGoil et al.ll20liTl . 
PRIMUS is a spectroscopic faint galaxy redshift survey 
to z ^ 1 using the Inamori Magellan Areal Gamera and 
Spectrograph camera on the Magellan I Baade 6.5 m 
telescope at the Las Campanas Observatory. PRIMUS 
observed ^ 2500 objects at once over a 0.18 deg^ field 
of view employing a low-dispersion prism and slitmasks. 
PRIMUS has covered a total of 9.1 deg^ of sky to a depth 
of iAB ^ 23.5 in seven well studied fields, one of which 
is the XMM-LSS field. It is in this field that the highest 
number of PRIMUS spectroscopic redshifts have been 
obtained. The redshift distribution of PRIMUS sam¬ 
ple galaxies peaks at z 0.6 and extends to z = 1.2 
for normal galaxies (the galaxies of interest in this pa¬ 
per) and z = 5 for broad-line active galactic nuclei. In 
the XMM-LSS field, PRIMUS observations cover ^ 2.88 
deg2 sky area and detected a total of 102218 objects. Due 
to the usual limitations of a prism based spectrograph, 
not all redshifts are equally robust. The PRIMUS team 
assigned a redshift confidence flag Q for every galaxy in 
their sample. They consider a redshift measurement to 
be robust if Q = 3 or 4. With this criterion, robust red¬ 
shifts were obtained for a total of 44451 normal galaxies 
and AGN in the XMM-LSS field. We only use the galax¬ 
ies with robust redshift measurements in our analysis. 
Note that in the XMM-LSS field there are also PRIMUS 
calib ration observations o f some sources in the VVDS 
field (iLe Fevre et al.l[2005ll . This was done to verify the 
quality of prism spectroscopy in PRIMUS by compar¬ 
ing it to higher resolution spectroscopy done with other 
telescopes. Since these calibration sources are mostly a 
mix of galaxies from different spectroscopic campaigns, 
it is difficult to perform any kind of statistical analysis 
on them. We have therefore excluded all sources in the 
calibration fields in the XMM-LSS area from our sample. 

2.1.2. Removal of AGN 

The PRIMUS color selection is designed to select nor¬ 
mal (non-AGN) galaxies with 0 < z < 1.2. Separately, a 
variety of t argeting criteria were used to target candidate 
AGN Isee iCoil et al.1 [201111 . The PRIMUS team care¬ 
fully identified the AG Ns in their sample by fitting AGN 
spectral templates Isee lGool et al.1 [201^ . These are in¬ 
dicated by the GLASS keyword in the PRIMUS catalog. 
We have removed the objects which have CLASS=AGN 
in our analysis. 

2.1.3. Removal of red galaxies 

A second, more serious, complication is that the 
PRIMUS sample contains both “blue cloud” and “red 
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TABLE 1 

Multi waveband surveys of the XMM-LSS field. 


Survey 

Total Area 
(deg2) 

Resolution 

X ") 

5(J sensitivity 
(mJy) 

A^tot 

A^primus 

-Nmatch 

GMRT 0.325 GHz 

12 

9.4 X 7.4 

0.75 

3929 

894 

111 

VLA 1.4 GHzi 

1.3 

5x4 

0.10 

505 

478 

109 

SWIRE 24Atm2 

10.6 

5.6 X 5.6 

0.45 

24799 

9231 

1812 

SWIRE 70Atm2 

10.4 

16.7 X 16.7 

2.75 

802 

301 

76 

SWIRE 160/im2 

10.3 

35.2 X 35.2 

17.5 

286 

106 

28 

HerMES 250/rm® 

18.87 

18 X 18 

25.8 

37905 

7331 

3004 

HerMES 350/rm® 

18.87 

25 X 25 

21.2 

42398 

8361 

3378 

HerMES 500/rm® 

18.87 

37 X 37 

30.8 

36933 

7293 

2856 


Atot: Total number of catalog sources - ^Simpson et alJ (|20Q6l L ^Surace et al[3; ^Smith et al.l l|2Q12ll : [Roseboom et al.1 l|201Qll 
-^PRIMUS' Number of sources within the PRIMUS footprint 
Amatch- Number of sources having a PRIMUS counterpart 

^http: //irsa.ipac.caltech.edu/data/SPITZER/SWIRE/docs/delivery_doc_r2_v2. pdf 



Fig. 1. — Footprints of the various surveys of the XMM-LSS field used in this work overlaid on the 0.325 GHz GMRT mosaic in units 
of Jy beam“^. We have saturated the image at 5 mjy bearn”^ for better representation. The grids have a separation of 1 degree in each 
axis. The thick black box out lines the 12 deg^ regi on having smooth rms noise of 150 /iJy beam“^. The red area shows the VLA 1.4 
GHz coverage of 1.3 deg^ from ISimpson et al.r(l200Q) . The inner blue area shows the 2.88 deg^ coverage of the PRIMUS and the outer blue 
area shows the coverage of CFHTLS. The yellow area outlines the SWIRE coverage of ~ 9 deg^ and the orange area shows the HerMES 
coverage of the XMM-LSS field covering ~ 15 deg^. 
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sequence” galaxies. These correspond to different galaxy 
populations with large differences in physical parame¬ 
ters such as the current star-formation rate, star forma¬ 
tion history, stellar mass, dust content and to a large ex¬ 
tent, galaxy morphology (jToieiro et al.ll201^ . Broadly, 
the “blue clo ud” galaxies represent active star-forming , 
disk galaxies (|Kauffmann et al.ll2003t iWvder et al.ll2007r) 
while the “red sequence” galaxies are p a.ssive and/or 
dust-obscured active sta r-forming galaxies (iBaldrv et al.l 
120061: iTavlor et al.ll201ljl . The significant and systematic 
differences in star formation properties of these two pop¬ 
ulations are likely to impact their radio and FIR prop¬ 
erties as well. It is therefore important to remove the 
“red” galaxies from our sample. 

Since the radio-FIR correlation is believed to be driven 
by star formation, we focus our attention on studying the 
“blue cloud” galaxies (henceforth referred as blue galax¬ 
ies). Blue galaxies are widely distributed across red shifts 
(see e.g., iLabbe et al.ll200'i^ iBrammer et al.l [201 Ih . ac¬ 
tively forming stars and are supp osed to dominate th e 
cosmic star formation rate density (|Magnelli et al.l[2009f) . 
Further, they are sufficient in number to produce bet¬ 
ter signal-to-noise ratio (S/N) for stacking analysis. The 
number of red galaxies in the PRIMUS sample are insuf¬ 
ficient for a robust analysis and are therefore not studied 
here. We find a total of 36,776 blue galaxies in PRIMUS 
at whose positions we perform image stacking in various 
infrared and radio bands. 

A color-magnitude diagram is used to separate 
“blue” and “red” galaxies. For the PRIMUS sample, 
iSkibba et al.l (1201411 proposed a color-magnitude diagram 
based separation using the u — g color and the g—band 
absolute magnitude {Mg). T he dividing line in th e color- 
magnitude plot is defined as (|Skibba et al.ll20ll l. 

{u — g)cut = —0.031Mg — 0.065z -I- 0.695. 

Here, z is the redshift of the galaxy. The redshift de¬ 
pendence accounts for the change in the positions of 
the two populations with redshift. Using this dividing 
line, we separate the PRIMUS galaxies into “blue cloud” 
and “red sequence” galaxies. Overall ~80 percent of 
the PRIMUS galaxies are blue, however, this fraction 
changes slightly with redshift. 

A small fraction of dusty star-forming galaxies may 
appear red due to reddening and thus can be classified 
as red galaxies owing to the color selection. Our sample 
misses star-forming galaxies where the blue light is sig¬ 
nificantly absorbed or sca ttered by dust. For PRIMUS 
galaxies. iZhu et aP (|2011ll showed that, at L ^ A*, ob¬ 
scured star-forming galaxies comprise ^15 percent of the 
red sequence population over all redshifts. They find, at 
lower luminosities (up to 0.2A*), the fraction of obscured 
star-forming galaxies in the red sequence population are 
up to ^ 30 percent. In the PRIMUS sample, due to the 
inherent flux limitation, such galaxies reside mostly in 
the lower redshift bins. At higher redshifts the low lu¬ 
minosity galaxies are missed by the PRIMUS itself and 
therefore it is difficult to study the overall fraction of 
obscured star-forming galaxies. However, in our sample 
of PRIMUS galaxies we find ~ 20 percent to be red, in 
the redshift range 0-1.2. Thus, up to 6 percent of the 
entire PRIMUS sample studied here could be dusty star¬ 
forming galaxies which are missed in our sample of blue 
galaxies. 



Fig. 2.— A small portion of our 0.325 GHz image is shown with 
the blue PRIMUS galaxies overplotted as red circles. It is apparent 
that the surface density of PRIMUS sources far exceeds the density 
of radio sources. Overall, 99.7 percent of the PRIMUS galaxies do 
not have detected 0.325 GHz radio counterparts above 5cr level. 

2.1.4. Malmquist bias 

PRIMUS is a flux limited survey. It samples 100 per¬ 
cent of target galaxies to iab = 22.5 and samples galaxies 
sparsely up to 0.5 magnitudes fainter with well defined 
a priori sampling rates. Like all flux limited surveys, it 
is subject to the well known Malmquist bias. At pro¬ 
gressively higher redshifts, galaxies with higher intrinsic 
luminosities form the luminosity cutoff in the sample. 
This bias could be corrected by employing a luminosity 
cutoff that is sufficiently high. Unfortunately, given our 
requirement for large sample size while stacking, such an 
approach would not work for us. For a fixed flux limit, 
and with no fc—correction, the ratio of the minimum lu¬ 
minosity in the highest redshift bin to that in the lowest 
one is about 125, for our chosen cosmological parameters. 
The lowest redshift bin samples galaxies that are more 
than ^ 0.2 times as luminous as the Milky Way while 
the highest bin samples galaxies that are more than ^ 25 
times brighter than the Milky Way. This implies that a 
true comparison of similar type of object (e.g. at a fixed 
fc—corrected luminosity or at fixed stellar mass) across 
redshifts is not possible for our sample. 

Further, PRIMUS is forced to exclude regions around 
bright stars, so there are gaps in the spectroscopic cover¬ 
age where such stars are present. Also, the spectrograph 
places limits on the minimum separation between slits. 
This implies that the highly clustered component of the 
sample is incompletely sampled. 

Despite these limitations and biases, we chose to work 
with the PRIMUS since it provides the largest sample of 
flux-limited, high-density, non-AGN galaxies with spec¬ 
troscopic redshifts over a relatively wide area to z ^ 1.2 
that is currently available. 

2.2. Data 

2.2.1. Radio data at 0.325 GHz and 1.4 GHz 
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The XMM-LSS field is observed at 0.325 GHz using 
the GMRT for 40 hours divided into four 10 hour obser¬ 
vation runs. The field is covered by a 16-pointing mosaic 
with a phase center to phase center separation of 1.0 de¬ 
gree. Scans are carried out in semi-snapshot mode of 6-17 
min each to optimize the itw-coverage. This resulted in 
approximately uniform sensitivity within the central re¬ 
gion. The final map covers the central ~ 12 deg^ having 
an average la rms noise of ~ 150 ^Jy beam“^ (shown 
as the thick black region in Figured]). Additional area 
covered has higher rms noise. A detailed description of 
the observations and data reduction will be presented in 
a forth-coming paper (Sirothia et ah, in prep). 

In our 0.325 GHz observations, we detected 3929 
sources with signal-to-noise ratio higher than 5 a local 
rms noise within the central region of which 894 sources 
lie within the PRIMUS footprint. There are only 111 
sources having a PRIMUS counterpart (see Table [T|) and 
85 of them are identified as blue galaxies. Thus, more 
than 99.7 percent of the PRIMUS galaxies are not di¬ 
rectly detected at 0.325 GHz. The large number of 
non-detections in the radio, makes this sample ideal for 
stacking analysis. Figure [2] shows the positions of the 
PRIMUS blue galaxies (red circles) overlaid on a small 
part of the 0.325 GHz radio image. The image clearly 
shows that bulk of the radio sky is blank at the position 
of the PRIMUS blue galaxies. 

To estimate the spectral index for fc—correcting the 
0.325 GHz luminosities to respective rest frequency at 
higher redsh ift we used the Very Large Array 1.4 GHz 
image from iSimpson et al.l (1200611 (shown as the red 
area in Figure [T]). They observed ~ 1.3 deg^ of the 
Subaru/XMM-Newton Deep Field and produced a cat¬ 
alog of 505 sources covering the central 0.8 deg^ with 
peak flux density limit of 100 ^Jy at 5cr. There are 478 
sources lying within the PRIMUS coverage of which 109 
have a PRIMUS counterpart. 

2.2.2. Fin data between 24 — 500 /rm 

The Level 6 Herschel Multi-tiered Extragalactic Sur- 
ve}{3 (HerMES) in the XMM-LSS field cover ~ 19 deg^ 
having 5cr sensitivity of about 25.8, 21.2 and 30.8 mJy 
at 250, 350 and 500 utur resp ectively using the SPIRE 
instrument (|01iver et al.ll201^ . The footprint of the Her¬ 
MES is shown in orange in Figure [T| HerMES provides 
images having angular resolution ~ 18, 25 and 37 arcsec 
at 2 50, 350 and 500 /rm, respectively. The p R2 cata¬ 
log ([Smith et al.l 1201^ iRoseboom et al.l 1201011 contains 
44120 sources in the XMM-LSS field that have been de¬ 
tected in at least one of the SPIRE bands of which 31676 
sources have been detected in all three bands. Table [T] 
lists the number of sources detected in each of the Her¬ 
MES bands along with the number of sources within the 
PRIMUS coverage. Only 2445 of the PRIMUS galaxies, 
i.e., about 5.5 percent, are detected in all the three bands 
of the HerMES. 

This field is also covered by the Spitzer space telescope 
as a part of Spitzer Wide-area InfraRed Extragalactic 

^ The HerMES project (http://hermes.sussex.ac.uk/) is a Her- 
schel Key Programme utilizing Guaranteed l ime trom the SPIRE 
instrument team, ESAC scientists and a mission scientist. The Her¬ 
MES data was accessed through the Herschel Database in Marseille 
(HeDaM - http://hedam.lam.fr) operated by CeSAM and hosted 
by the Laboratoire d'Astrophysique de Marseille. 


(SWIRE) survejl3 (|Lonsdale et al.|[2003l using IRAC 3.6, 

4.5, 5.8 and 8 fim bands and MIPS 24, 70 and 160 fim 
bands over ~ 10 deg^ (see Figure [T}. For modeling the 
re-processed dust emission in the infrared bands, we com¬ 
bined imaging from the SWIRE0 at A = 24, 70 and 160 
^m with the HerMES DR2 at 250, 350 and 500 p,m. 

3. ANALYSIS 
3.1. Image stacking 

Most radio images are overwhelmingly dominated by 
“blank sky” because radio source densities on the sky are 
very low, at the sensitivity we can achieve with current 
technology. Even in our relatively deep 0.325 GHz image 
we only detect about 500 radio sources per square degree 
above ^ 500 /rJy (> 3cr). In contrast, the dense sampling 
in PRIMUS yields ^ 2500 galaxies per 0.18 deg^, corre¬ 
sponding to a source density of ^ 14000 deg“^. In such a 
situation, stacking on the radio image at the positions of 
the PRIM US galaxies allows us to extract signals from 
the noise (|White et al.l 12007]) . The stacking process is 
straightforward - it involves making a cutout centered 
at the position of each sample galaxy and then obtaining 
a mean stack image by averaging over all cutouts on a 
pixel-by-pixel basis. 

We note that some authors have chosen to do a median 
stack instead of a mean (average) stack le.g. lWhite et al.l 
l2007f) . The median has the advantage of not being bi¬ 
ased by a few outliers anywhere in the image. However, 
the interpretation of the median value for low S/N data 
typically involved in image stacking is complicated. For 
such data, simulations carried out bv I White et al.l (120071) 
show that the computed median value is shifted from the 
true median towards the “local mean” value. Here, the 
“local mean” is defined as the mean of the values within 
approximately one rms of the median. Unfortunately, 
the degree of the shift depends on the noise level; as the 
noise increases, the recovered value approaches the local 
mean. Due to this effect, the recovered median value is 
a function of both the intrinsic distribution of the pa¬ 
rameter bei ng measured ( in our case, radio flux) and the 
noise level. iWhite et al.l (1200711 note this limitation but 
choose to go with the median stack since their sample - 
a large number of quasars drawn from the Sloan digital 
sky survey - has quite a few extreme outliers in the form 
of radio-loud quasars. Our sample is free of extremely 
bright radio counterparts; in fact, more than 99 percent 
of our sample galaxies have no detectable radio counter¬ 
part at 0.325 GHz even with our deep radio imaging. For 
this reason, we chose to perform mean stacks since the 
statistical properties of the mean are very well under¬ 
stood, even in low S/N situations. 

We stacked the images at 0.325 and 1.4 GHz in the 
radio and at 24, 70, 160, 250, 350 and 500 ^m in the in¬ 
frared at the positions of the blue galaxies. We divided 
the sample into redshift bins of 0.1 between redshifts 0 to 
1.2. Table O lists the number of galaxies in each redshift 

^ The Spitzer Space Telescope is operated by the Jet Propul¬ 
sion Laboratory, California Institute of Technology, under contract 
with NASA. SWIRE was supported by NASA through the Spitzer 
Legacy Program under contract 1407 with the Jet Propulsion Lab¬ 
oratory. 

® Downloaded from: http://swire.ipac.caltech.edu/swire/astronomers/ 
data_access. html 




























6 


A. Basu et al. 


bin for each of the datasets used in this work. Stack¬ 
ing was performed separately in each redshift bin. The 
PRIMUS sample has large number of galaxies in each 
redshift bin with a minimum of 1100 galaxies in the bin 
1.1 to 1.2 and a maximum of ^ 5300 galaxies in the bin 
0.6 to 0.7. The large number of galaxies in each bin al¬ 
lows us to detect signal from extremely faint galaxies in 
the stacked images with high signal-to-noise ratio (> 4). 
Due to smaller sky coverage at 1.4 GHz, the number of 
PRIMUS blue galaxies available for stacking on the 1.4 
GHz image in each redshift bin, are lower than those at 
other wavebands. 

Before stacking, a 61 x 61 pixel-image centered at the 
position of each galaxy within the redshift range of the 
corresponding bin was cutout from the images of the sur¬ 
veys mentioned in Section 2. The flux density image of 
each cutout was converted into a luminosity image using 
the respective redshift before doing a mean stack. We 
used the peak luminosity in the stacked image at the po¬ 
sition of the PRIMUS sources (typically the center pixel) 
to e stimate the mean lu minosity at that redshift bin (see 
e.g. iWhite et al.l I2007ll . For unresolved point sources, 
the peak emission is the same as the total emission in¬ 
tegrated within the resolution beam. However, this may 
differ i f the sources being stacked are extended w.r.t the 
beam (|Karim et al.ll201l]l . In our study, the PRIMUS 
galaxies are point sources for ^ 9 arcsec beam at 0.325 
GHz. Hence, the peak-pixel value of the stacked imaged 
can be used as the integrated luminosity of the sample. 

We carried out stacking of the flux density images 
also. The peak value of the luminosity or flux den¬ 
sity in stacked images should be computed w.r.t zero 
mean background. However, due to the contribution of 
emission from other sources, the background level in the 
stacks is not zero mean and is slightly offset towards posi¬ 
tive values. We have therefore subtracted this offset value 
from all pixels in the stacked images at each redshift bin. 
The mean flux density, the stacked image rms noise and 
the mean luminosity for each redshift bin are tabulated 
in Table El 

We note that the luminosity derived from the stacked 
fluxes and mean redshifts differ from the stacked lumi¬ 
nosities especially in the low redshift bins (see Table El). 
We list the luminosities obtained from stacking in lu¬ 
minosity space (column 7) and from stacked fluxes and 
mean redshift (column 8). In the lower redshift bins 
{z < 0.3) they differ significantly by up to ^^30 percent. 
However, at higher redshifts {z > 0.3) the two meth¬ 
ods yield similar values within 5 percent. This is mainly 
because the luminosity at the edges of the lowest 2 bin 
(0.1 < 2 < 0.2) varies by a factor of 4.5 as compared to 
a factor of 1.6 for a higher z bin (1 < z < 1.1). The 
stacked fluxes in a redshift bin do not account for this 
variation and therefore we chose to perform stacking in 
luminosity space. 

For robust error estimation of stacked luminosities 
and fluxes, we performed standard bootstrap analysis 
(|Efron fc Tibshiranilll994l ) on the cutout images in each 
redshift bin. At each redshift bin, a cutout was ran¬ 
domly selected A^obj times (allowing repetitions) and 
were stacked. Here, Ahbj is the number of objects in 
a given redshift bin (see Table El- We performed 20 
bootstrap iterations and computed the mean luminosity 
and the flux density for each. The standard deviation 
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Fig. 3.— Radio spectral index between 0.325 GHz and 1.4 GHz 
with redshift calculated using the mean stacked luminosities. The 
dashed line shows the typical spectral index of 0.8 observed for 
nearby normal galaxies. The errors are at Icr. 

for these 20 iterations is considered as the error on the 
stacked quantities at a given redshift bin. 

Determination of the luminosity of the stacked images 
could be a ffected if the sources in th e primary sample are 
clustered (IBethermin et al.l 1201011 . ISkibba et al.l (|2014fl 
found evidence of clustering of the PRIMUS galaxies, 
especially for the red galaxies, which can affect the es¬ 
timated stacked luminosity for images with coarse reso¬ 
lution , like that of the HerMES bands. IBethermin et al.l 
(|2012t) showed that the contribution due to clustering can 
be estimated by modeling the radial profile of the stacked 
image as a point spread function {psf) and convolution of 
the psf and two-point angular correlation function of the 
parent sample. In our case, at all the wavebands, the 
radial profile can be fitted by a single component Gaus¬ 
sian profile within 7 percent, that represents the psf even 
at 500/rm where the resolution is poorest. This indicates 
that the correction due to clustering of the PRIMUS blue 
galaxies is less than 10 percent and is within the errors 
estimated using the bootstrap method, and hence does 
not significantly affect our results. Thus, working with 
blue galaxies has the advantage of being significantly less 
clustered and such corrections are small enough to be ig¬ 
nored. 

3.2. Spectral energy distribution and k—correction 

We obtain rest-frame emission by applying 
fc—correction to the observed luminosities in radio 
and FIR bands. In this study, we explore the radio-FIR 
correlation between rest-frame radio emission at 1.4 
GHz and FIR emission both monochromatic at 70/rm 
and bolometric between 8-1000/rm. These specific radio 
and FIR bands are chosen to en able direct comparison 
with the literat u re (see e.g. . lAppleton et alJ l2004 
Mao et all 120111: H vison_et ahl I2010allb| : iBourne et al l 
20111 IMagnelli et al.11201^ In the radio, estimating 
the luminosity at a given rest frequency, in our case 
1.4 GHz, is relatively easier owing to power law nature 
of the spectrum. Figure El shows the variation of the 
spectral index, Ont, defined as Si, oc determined 

using the mean stacked luminosity at each redshift bin 
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Fig. 4.— The fitted SED of the stacked sources in infrared bands 
at various redshift bins shown in colors. The SEDs are fitted jointly 
with a mid-IR power law and a single temperature greybody. 


between 0.325 GHz and 1.4 GHz. The spectral index 
in each redshift bin is consistent within Itr errors with 
typical spectral index of ^ 0.8 (shown as dashed line 
in Figure |3l) observed for n e arby normal star-f orming 
galaxies (|Niklas et al.l Il997t iBasu et ^ l2012all . We 
therefore estimated the rest frame 1.4 GHz luminosity 
(l'1.4GHz,rest) by extrapolating the observed frame lumi¬ 
nosity at 0.325 GHz (To. 325 GHz,obs) assuming the typical 
non-thermal spectral index of 0.8 for normal galaxies. 

Thus, Tl. 4GHz,rest = ro. 325 GHz,obs [0.325(1 -|-z)/1.4]° ®, 
where z is the redshift. 

At FIR wavebands, due to the complex nature of the 
spectral energy distribution (SED), robust modeling is 
necessary for accurate estimation of the monoc hromatic 
and bolometric FIR luminosities. iCasevI (|20I2f) pointed 
out that htting a simple modified blackbody spectrum 
(henceforth referred to as greybody) to the FIR spec¬ 
trum is inadequate in fitting the mid-IR observation, 
whereas htting a multi temperature greybody spectrum 
introduces several free parameters. Further, they demon¬ 
strated that a model consisting of single temperature 
greybody + midHR pow er law hts the FIR spectrum well. 
Following iGasevI (1201211 we modelled the FIR spectrum, 
S{X), as. 


S'(A) = Agb 


(1 -e-^^)A-3 

^^hc/XkT _ 2 ^^ 


+ ApL 


= -(A/Ac)" 


( 1 ) 


Here, Aqb and App are the greybody and mid-IR power 
law amplitude normalization respectively, Ac is the mid- 
IR turnover wavelength, a is the mid-IR power law in¬ 
dex, T\ = (Ao/A)^ is the optical depth and is unity at 
Ao (assumed to be 200 ^m), (3 is the dust emissivity in¬ 
dex, T is the greybody temperature and h, c, k are the 
Planck constant, speed of light and Boltzmann constant, 
respectively. Characteristic dust temperature (Tdust) is 
given by Wien’s displacement law, Tdust = ^7 Apeak (Mm), 
where, b = 2.898 x lO^^m K and Apeak is the peak wave¬ 
length of the htted SED. 


Although there are six htting parameters, namely, 
Agb, ApL, a, /?, T and Ac, the parameters App and 
Ac are coupled to the rest, ther eby reducing the number 
of free parameters to four Isee iCas^ 12012 1. Thus, for 
a robust htting of the EIR SED, at least 5 data points 
are necessary. We have therefore stacked, in luminos¬ 
ity space, 6 wavelength bands in the EIR, i.e., three 
MIPS bands (24, 70 and 160 ^m) from SWIRE and three 
SPIRE bands (250, 350 and 500 ^m) from HerMES at 
the position of the PRIMUS galaxies to constrain the 
SED. However, to constrain the mid-IR power law index, 
a, at least 3 photometric points in the mid-IR (between 
^10 — 50 /rm, rest-wavelength) are required, which are 
not available in our case. We have therefore hxed a = 2 
Sirn i larly, we hx ed (3 = 1.5 (following 
iMagnelli et al.l[2?)T^ ICasevI[201^ . 

Separate hts were done for the stacked SEDs in lumi¬ 
nosity space, in each of the 12 redshift bins shown in Eig- 
ureSl Each ht is color coded for the redshift bin. In Eig- 
ure[8]we show the two components, greybody (red dashed 
curve) and mid-IR power law (blue dash-dot curve) at 
each redshift bin along with reduced chi square (Xred) 
and errors in the htted parameters. The two compo¬ 
nents together give an excellent ht to the data with one 
exception. The SED ht for the hrst redshift bin is poor 
having large and more than 50 percent error in the 
derived EIR luminosity. We have excluded this bin from 
all our further analysis. 

To compare our results using stacking, and study its 
nature in perspective of the large span in luminosity (^ 5 
orders of magnitude), we also study the radio-EIR cor¬ 
relation using detected sources along with the stacked 
sources. Eor the sources that are detected at 0.325 GHz, 
we cross matched them with the HerMES DR2 catalog 
(iRoseboom et al.l[20I0ll . There are 1805 sources that are 
detected at 0.325 GHz and in all the three HerMES bands 
(250, 350 and 500Mm). The HerMES bands are required 
for a better constraint on the “cold dust” greybody spec¬ 
trum. However, we note that for the above SED fitting 
method to work best in estimating the TO^m monochro¬ 
matic luminosity, one needs mid-IR observation to con¬ 
strain the power law component. We therefore, addition¬ 
ally imposed the condition, that the sources be detected 
at 24 and 70 from the SWIRE survey. There are 
231 sources detected in hve infrared bands from 24^m to 
bOO^m with 0.325 GHz counterparts. Only 26 of the 231 
sources have spectroscopic redshift from the PRIMUS, 
predominantly due to small areal coverage. Therefore, 
to increase the sample size we used photometric red- 
shifts (zphot) from the Ganada-France-Hawaii Telescope 
Legacy Surve}!! (CFHTLS; lllbert etldl 120061 1 available 
in the XMM-LSS field that covers a larger area (see 
Figure [T|). GFHTLS provides accurate Zphot for galax¬ 
ies, with only ~ 4 pe rcent catastrophic failures, i.e., 
Az/(1 -I- Zspec) > 0.15 (jllbert et al.ll2006[l . Of the 231 

^ The CFHTLS public data release includes observations ob- 
tained with MegaPrime/MegaCam, a joint project of CFHT and 
CEA/DAPNIA, at the Canada-France-Hawaii Telescope (CFHT) 
which is operated by the National Research Council (NRC) of 
Canada, the Institut National des Science de TUnivers of the Cen¬ 
tre National de la Recherche Scientifique (CNRS) of France, and 
the University of Hawaii. The data products are produced at TER- 
APIX and the Canadian Astronomy Data Centre, a collaborative 
project of NRC and CNRS. 
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Fig. 5.— Left hand panel: Luminosity at 1.4 GHz (Li, 4 GHz) versus monochromatic luminosity at 70/rm (Lyo/im) at rest-frames. The 
stacked sources are shown as stars and the sources detected in the XMM-LSS field are shown as circles. The symbols are color coded based 
on their redshift. The solid line shows the fit to the entire data, the dashed and dashed-dot lines are for the stacked and detected sources 
respectively. The slopes are listed in Table 2. Typical errors on the data are shown in the lower right. Right hand panel: Li. 4 GHz versus 
bolometric luminosity (Ltir)- The symbols and the lines have the same meaning as in the left hand panel. 


sources, 126 have 2;phot available. These 126 sources are 
an assortment of galaxies comprising of 100 blue and 26 
red galaxies spanning up to a redshift of 0.95. Note that, 
19 of these 100 blue galaxies lie in the PRIMUS coverage 
and have also been included while stacking. 

The choice of a plays a crucial role in determining the 
ApL and Ac and, in turn the SED fit below 100 /tm, 
leading to deviant fits. We have therefore performed a 
grid search for the value of a between 1.2 and 3.2 in steps 
of 0.05. We fixed a where the to the fit was minimum. 
More than 90 percent of the objects could be fitted with 
a in the range 1.8-2.3. 

We would like to emphasize that, the three HerMES 
data points help in constraining the “cold dust” greybody 
spectrum while for the mid-lR power law part, more than 
one data point below ^ 100 /im observed-frame wave¬ 
length (in our case 24 and 70 ftm) is necessary for robust 
estimation of fc—corrected 70/rm luminosity. Using only 
24/im and three HerMES bands, constrains the FIR grey¬ 
body well, but the mid-lR SED is not well constrained 
leading to increased dispersion in the radio-FlR correla¬ 
tion when studied at 70/im. Having the observed frame 
70/rm data point is crucial to our analysis. 

4. RESULTS 

Image stacking technique allows us to estimate mean 
luminosities of our sample galaxies in each redshift bin 
of 0.1 across 0 < z < 1.2 (see Section 3). We 
summarize our results of the stacked images at var¬ 
ious wavebands in Table [31 Note that, we detect 
emission from faint galaxies that have Ti, 4 GHz of few 
times ~ 10^° W Hz“^ and of few times 10^^ 

W Hz“^. Using stacking we probe radio-FIR correla¬ 


tion for galaxies that are up to 2 order to magnitude 
fainter than the samples used in previous studies based 
on direct detections in the similar redshift range (see e.g . 
Garretd [2002t [Appleton et al.l I2004 iKovacs et al.l 120061 : 


Mao et al.l l2011l : lDel Moro et al.ll2013ll . However, this 
increased depth is restricted to the lowest redshift bins, 
due to the Malmquist bias present in the parent optical 
sample. The radio/FlR stacking process does not intro¬ 
duce any additional biases, but it cannot overcome the 
inherent bias in the parent optical sample. 

The radio-FIR correlation is generally quantified by 
two parameters - 1) the slope in log-log space, b, given 
by, Tradio oc and 2) the parameter ‘g’ dehned 

as q = logio(AiR/Aradio)- Here, Aradio is the lumi¬ 
nosity at a radio frequency, in our case at 1.4 GHz 
rest-frequency. However, several definitions of Air are 
found in the literature, such as rest-frame monochro¬ 
matic lu minosity at some specified infrared w avelength 
(see e.g. lAppleton et ah! 12004 iMao et aL[l2011h or bolo¬ 
metric l uminosity integrat ed between 40 — 120 ftm. (de- 
fine d by iHelou et al.||1985 D or betw e en 8 — 1000 iim (see 
e.g iBeill 20031: Ivison et aP I2010bl : iBourne et al.l 120111 : 
iMagnelli et al.l[2014f) . Here we study the radio-FIR cor¬ 
relation using both rest-frame monochromatic infrared 
luminosity at 70/rm (Ayo^m) as well as bolometric to¬ 
tal infrared luminosity integrated between 8 — 1000 /rm 
(Atir)- For the monochromatic case, we define qropm = 
logio(A 70 ijm/Ai, 4 GHz) and for the bolometric case, we de- 
Hne, 


?TIR = logio 


Atir(W) 
3.75 X 1012 


logio [Al. 4GHz (w Hz ^)] 

( 2 ) 
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Fig. 6.— Left hand panel: The variation of qro^m with redshift (1 + z). The stars and circles represent stack ed and detected source s 
respectively. The solid and the dashed lines shows the mean value of qro^m and Icr dispersion respectively from [Appleton et al.l 1120041) . 
Right hand panel: The variation of qfir with redshift (1 + z). The red dashed line shows the fitted curve for the stacked sources of the 
form (2.53 ib 0.04)(1 + The solid line shows the mean value of Qtir from lBell (120031) observed for local galaxies. The mean 

values of these quantities obtained in our study are given in Table 2. 


follow ing iHelou et al.l (jl985[ l. Note that iHelou et alJ 
(119851) defined the bolomet ric infrared luminos ity inte¬ 
grated between 40 — 120 /rm. iBourne et alJ (I2011D showed 
that the use of 8 — 1000 luminosities increases the 
absolute value of qTm by ^ 0.32 corn pared to that orig¬ 
inally defined bv IHelou et alJ (jl985f ). The rest-frame 
far-infrared monochromatic luminosity at 70/xm and the 
bolometric luminosity are computed using the SED fits 
discussed in Section 3.2. 

In Figure [5] (left hand panel) we plot the radio lumi¬ 
nosity at rest frame 1.4 GHz (ii. 4 GHz) vs. the FIR lumi¬ 
nosity at 70/rm (TTO/^m) for the stacked galaxies (shown 
as stars) and the directly detected galaxies (shown as cir¬ 
cles). The symbols are color coded based on the redshift. 
We note that the mean stacked luminosities in the higher 
redshift bins are higher. This trend is expected due to the 
inherent Malmquist bias. The radio and FIR luminosities 
are found to be strongly correlated with Spearman’s rank 
correlation, r > 0.99. We fitted for the slope of the cor- 
relation using ord inary least-square “bisector method” 
(|lsobe et al.lll99Clll in log-log plane. The solid line shows 
the fit to all the data points and have a slope 1.09±0.05. 
We also did separate fits for the stacked and detected 
galaxies. The dashed line shows the fit for the stacked 
galaxies for which the slope is found to be 1.04±0.03. 
The dash-dot line represents the fit for detected galaxies 
and have a slope of 1.12±0.05. They all agree within er¬ 
rors suggesting the radio-FIR correlation holds over the 
entire luminosity range of five orders of magnitude, prob¬ 
ing galaxies with different ISM properties. The values of 
the slopes are tabulated in Table 2. 

Figure [5] (right hand panel) shows the radio-FIR cor¬ 
relation between Ti, 4 ghz and the total infrared luminos¬ 


ity (Ttir)- We find strong correlation between the two 
quantities. In this case, we find the slope to be signif¬ 
icantly steeper than unity. The slope of the radio-FIR 
correlation is found to be l.II ± 0.04 when fitted for 
the stacked sources and detected sources together (solid 
line). By separately fitting the stacked galaxies and the 
detected sources we find the slopes to be 1.12 ±0.03 and 
1.13±0.04, respectively (see Table 2). Si milar non-lin ear 
slope of 1.10 ± 0.04 was also reported bv I Belli (I2003D for 
normal star-forming galaxies in the local universe having 
Ltir in the range 10® to 10^^ Lq. 

In Figure [HI we study the variation of q 70 fj,m (left hand 
panel) and qxiR (right hand panel) with redshift (1 ± z). 
The mean values of the parameter are tabulated in Ta¬ 
ble 2. We find the mean value of q 70 /jm to be 2.18 with Icr 
dispersion of 0.26. This is consistent with (qro^m) = 2.15 
(shown as the solid line) and Icr di spersion of 0.15 ( shown 
as the dashed lines) observed by lAppleton et al.l (1200411 
over similar redshift range and bv IMao et al.l (|201ll) for 
redshifts up to ^ 3 for more luminous galaxies. How¬ 
ever, for the stacked so urces, qroiim is f o und t o be lower 
than that observed bv lAppleton et all (|2004ll and have 
(qroiim) of 2.0 ± 0.2. Within errors, qroiim is observed to 
remain constant with redshift. 

qtir shows a mean value of 2.50 with a Icr dispersion 
of 0.24. This is clos e to the value 2.64 observed for local 
galaxies (|Bellll200^ with similar dispersion. The mean 
value of qTiR for the stacked sources are found to be 
lower, (^tir) = 2.34 ± 0.22, than that of the detected 
sources, (^tir) = 2.51 ± 0.24 (see Table 2). Our esti¬ 
mated value of ((/tir) for stack ed sources is close to the 
value 2.40 ± 0.29 observed by llvison et al.l ()2ni0aD for 
sub-mJy radio emitting galaxies. Note that, the radio- 
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TABLE 2 

Results of the radio-FIR correlation. 


70/jm TIR FIR 



Slope 

(grOMm) 

Slope 

(gTin) 

Slope 

(9fir) 

Stacks 

1.04±0.03 

2.0±0.3 

1.12±0.03 

2.34±0.22 

1.07±0.03 

2.07±0.21 

Detections (blue) 

1.05±0.05 

2.22±0.27 

1.12±0.06 

2.51±0.24 

1.07±0.06 

2.25±0.24 

Detections (red) 

1.30±0.13 

2.16±0.40 

1.14±0.07 

2.48±0.32 

1.29±0.12 

2.19±0.31 


Detections (blue+red) 1.12 ±0.05 2.22 ± 0.24 1.13 ±0.04 2.51 ±0.24 1.13 ±0.05 2.25 ± 0.24 

Stacks±detections 1.09 ± 0.05 2.18 ±0.26 1.11 ±0.04 2.50 ±0.24 1.10 ±0.04 2.23 ± 0.25 

The slopes are estimated by fitting a straight line in the log-log space of the luminosity plots shown in Figure [5] 
(...) represents the mean values of the quantity along with their measurement errors from Figure [6l 


FIR correlation is tighter when studied for bolometric 
FIR luminosity than for monochromatic 70/im luminos¬ 
ity. This is evident from the mean value of gxiR showing 
lesser spread around the mean (0.24 dex) as compared 
to 970 Aim (0.26 dex). In linear space this corresponds to 
^ 16 percent spread in Ti, 4 ghz vs. Ttir as compared to 
28 percent spread when studied between Li. 4 ghz vs. 

^70fim ■ 

For completeness with the various definitions of the pa¬ 
rameter ‘q\ we also present the results of the radio-FIR 
correlation between rest-frame Ti, 4 GHz and Lrir (listed 
in Table 2). Here, Trir is the bolometric infrared lu¬ 
minosity integrated between 40 — 120 ^m. There is no 
significant differences in the values of qpm and the slope 
when compared with that of monochromatic 70^m. This 
is not s urprising bec a use, as per the traditional defi¬ 
nition of lHelou et al.l (|1985ll . the factor 3.75 x 10^^ W 
(in Equation [2]) normalizes the integrated luminosity be¬ 
tween 40 — 120 /iiTL to the mean luminosity at ~ 80 /rm. 
This value is close to 70/Lm. 

4.1. Detected red galaxies 

To assess the effect of including the 26 red galaxies on 
the radio-FIR correlation, we have separately studied the 
red and blue galaxies that are directly detected. The two 
populations are well mixed in the radio-FIR correlation 
and there are no systematic difference between them. In 
Table 2 we list the slope and the ‘g’ parameter of the 
two classes estimated for TOfim, TIR and FIR. Both blue 
and red galaxies follow the correlation with similar ‘g’ 
parameters for all the three cases (i.e., gro/im, 9 tir and 
^fir) ■ However, the slope of the correlation is found to be 
steeper in red galaxies as compared to the blue galaxies 
when studied between Ti. 4 ghz in the radio and Lroiim 
and Tfir in the infrared (see Table 2). This difference is 
at less than 2a significance and is perhaps caused due to 
relatively smaller sample size of the red galaxies. We do 
not observe any difference in the slope for blue and red 
galaxies when studied for Ttir- Overall, including the 
red galaxies in the sample of detected galaxies does not 
bias our results. 

Note that in our stacking analysis we have not distin¬ 
guished between the objects that have been detected as 
a PRIMUS counterpart or not. This is crucial in esti¬ 
mating the typical mean flux or luminosity of the sample 
in each redshift bin and waveband. Since, we have less 
than 50 percent detections w.r.t PRIMUS objects, the 
median stacking could be insensitive to the objects hav¬ 
ing PRIMUS counterparts. In our analysis, we find that 
the mean quantities are at an average 20 percent higher 
than that of the median quantities. Our method of stack¬ 
ing is equivalent to the quantity measured in Equation 3 


of iMagnelli et al.l (|2014D . The difference between mean 
and median stacking is well captured in our estimated 
bootstrap errors. The use of median luminosities does 
not significantly affect the slope of the radio-FIR corre¬ 
lation shown in Figure [5] (left hand panel). However, it 
can systematically increase the value of qroam and gxiR 
of the stacks (the stars shown on Figure [ 6 |) by ~ 0 . 1 , 
i.e., ~ 5 percent. This change in and gxiR is well 

within the errors. 

5. DISCUSSION 

We have studied the radio-FIR correlation of an opti¬ 
cal flux-limited sample of “blue cloud” galaxies (Ttir ~ 
10 ® — lO^^L©) in the redshift range 0.1 — 1 . 2 , using image 
stacking technique in luminosity space, and for compar¬ 
ison, for luminous galaxies (Ttir ^ 1 O^^L 0 ) detected in 
the XMM-LSS Held using deep GMRT observations at 
0.325 GHz. IBourne et all (IWl . using a coarser redshift 
bin-size of > 0.2 for stacking, probed a similar luminosity 
regime for stellar-mass-selected galaxies. We measure the 
mean luminosity of blue galaxies to be at least 2 orders 
of magnitude fainter than the galaxies that are directly 
detected in radio and infrared in the XMM-LSS field and 
about three orders of mag nitude fainter than that of the 
previous studies (see e.g. , GarrettI 120021: 1 Appleton et al.l 
12004 iKovacs et al.l 1200^ Mao et al.l 1201111 . Our study 
probes a relatively poorly observed regime of the radio- 
FIR correlation in luminosity and in redshift. 

5.1. Variation of dust temperature 

In Figure [7] (left-hand panel), we show the variation of 
Tdust with the total infrared luminosity (Ttir in units of 
Lq) for the stacked sources (shown as blue stars) and the 
detected sources (shown as black points). Our result is 
consistent with the second order polynomial fit (shown 
as red dashed line) for slig htly more luminous g alaxies 
between 0.2 < z < 0.5 f r om IMagnelli et al.l (|2014li . Note 
that the IMagnelli et all (12014) curve does not extend to 
the faintest data point in our sample. 

In Figure [7] (right-hand panel) we plot the variation of 
the characteristic dust temperature, Tdust with redshift, 
1 + z. We find T^-ast to vary between 18 — 28 K with Tdust 
increasing linearly with redshift as, 

Tdust = (8.2 ± 0.9)(1 + z) + (11.5 ± 1.5) (3) 

A similar linear variation of T dust with redshift was re¬ 
ported bv iKovacs et al.l (|2006D for distant sub-millimeter 
selected galaxies which are > 4 orders of magnitude more 
luminous than our galaxies. However, such a trend could 
be caused due to Tdust being correlated with Ttir which 
in turn is correlated with redshift. 
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Fig. 7.— Left-hand panel Variation of dust temperature with the total infrared luminosity (Ltir) computed by integrating the SEDs 
between 8—1000 /rm. The blue stars represents stacked sourc es and circular points are for detected sources. The red dashed line shows 
the second order polynomial fit between 0.2 < z < 0.5 from iMagnelli et ahl II2014I1 . Right hand panel: Variation of the estimated dust 
temperature (Tuust) with (1 + z). The dashed line shows the linear best fit of the form Tdust = (8.2 ± 0.9)(1 + z) + (11.5 ± 1.5). 


5.2. Variation of ‘q’ 

The mean value of gxiR for the stacked sources is 
found to decrease with redshift, a nd can be modeled 
as ffTTR. PC (1 + z)'^ (see, e.g., iBourne et al.l 120111 : 
llvison et al.ll 2010 ^ . where 7 = —0.16 ± 0.03 (Figure [5] 
right-hand panel). The fitted value o f 7 is similar to 
7 = —0.18±0.10 found bv iBourne et al.l l) 20 111) for galax¬ 
ies with mass-lim it log(M ^) > 10.5 and 7 = —0.15 ±0.03 
by llvison et al.l (| 2010 al) for slightly luminous galaxies 
(±i. 4 GHz 10^^ — 10^® W Hz“^). This slight decrease 
in ( 7 tir with redshift, however, is not very significant 
within the large errors. The correlation coefficient was 
estimated taking into account the errors of gxiRi using 
a Monte Carlo method wherein 1000 random samples of 
<?TiR were generated within the Icr error for each redshift 
bin. The values of qtir were then randomly sampled 
from each redshift bin 1000 times to mimic different re¬ 
alizations and the Spearman’s rank correlation was es¬ 
timated for each. We found the mean rank correlation 
coefficient to be ~ —0.13. 

We do not observe such a trend for qjopm- This could 
be caused due to the rest-frame luminosity at 70/im be¬ 
ing dust temperature dependent. Figure [7] (right hand 
panel) shows an increase in dust temperature with red¬ 
shift resulting in increasing T 70 /im, hence compensating 
for any decrease in qro um- This trend is co nsistent with 
recent observations of ISmith et al.l (|2014D , where they 
find monochromatic q near the peak of the dust emission 
to remain constant with dust temperature. 

5.3. Is ‘q’ an indicator of evolution of the correlation? 

In our study, we find the radio-FIR correlation to hold 
true within the redshift range 0 . 1 - 1.2 for blue galaxies 
and also for much luminous galaxies spanning about 5 or¬ 


ders of magnitude. The correlation can be fitted by a sin¬ 
gle slope, b = 1.09 ± 0.05, when studied with monochro¬ 
matic FIR luminosity (T 70 /tm)- For bolometric FIR lu¬ 
minosity (Ltir), we find, b = 1.11 ± 0.04 (see Table 2 ). 
For our sample, we find that the slope of the radio-FIR 
correlation is systematically steeper than unity with l.Scr 
significance for Tro/im and 2.75a for Ttir- 
The non-linear slope of the radio-FIR correlation 
brings to light the ambiguity in use of the quantity ‘g’ 
(both g 70 /im and gxin) in quantifying the tightness of the 
correlation or in studying evolution of the radio-FIR cor¬ 
relation. Note that, apart from the traditional definition, 
q can also be written as, 

9 = - logio a + Ti.4GHz (4) 

Here, a is the proportionality constant of the radio-FIR 
correlation, i.e., Ti, 4 GHz = This relation holds for 

both monochromatic qro^im and bolometric gxiR- Clearly, 
for non-unity slopes, the absolute value of q would de¬ 
pend on the slope and also on radio luminosity and 
hence cannot be assumed to be a constant. Several 
other works i n the literature have indeed reported non¬ 
line ar slop es (iPrice &: Durid [TOOl: iNiklas &: BeddUO^ 
iBelillMil . Moreover, Ti, 4 ghz in the rest-frame depends 
on the assumed value of the spectral index and thus can 
affect the value of q further. It is, therefore, difficult 
to compare the value of the parameter q in the litera¬ 
ture and interpret any variation of q as evolution of the 
correlation. This non-linearity of the slope could also 
give rise to the observed slight change in gxiR (see Sec¬ 
tion 4 and Figure [ 6 |). We have modeled the variation as 
gxiR = 2.53 (1±2:)“°'^®, which results in Agxm ~ —0.26 
between z = 0 and 1. This change has been observed 
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in oth e r studies, e.g.. Ilyison et al.l (|2010ali ; iBourne et alJ 
(|2011f) : [Magnelli et alT i 20121) . Following Equation^ the 
expected change in gxiR is given by, 


A^tir = 


1 - b 


log 


10 


Li 

Li 


4GHz |z=l 
4GHz |z=0 


(5) 


Here, Ai. 4 GHz \z is the rest-frame 1.4 GHz luminosity at a 
redshift 2 . From Figure[5]we find Fi. 4 GHz to increase by a 
factor ^ 10 ^ between z = 0 and 1 for the stacked galaxies. 
Thus, AgxiR ~ —0.27 for b = 1.1 (see Table 2), similar to 
the observed change. It is therefore difficult to interpret 
any change in the parameter ‘g’ as an evolution of the 
radio-FIR correlation for non-linear slopes. Note that, 
although our results are based on a sample that is flux 
limited and color selected, it is important to ascertain 
the slope for a more representative sample to conclude 
any evolution based on ‘q’. 

The parameter ‘q’, explicitly depends on the absolute 
values of the ISM parameters, such as, density of dust 
(pdust), Tdust, dust emissivity (Q), properties of dust 
grains, like size (a), CRE density (riGRe), magnetic field 
(B) and the non-thermal spectral index (ant)- In terms 
of these quantities q can be written as; 

PdustQ(\a)B},(Td 

USt) /^\ 

«- 

Hence, based on a controlled sample ‘ 9 ’ can be used as 
a proxy to ‘relatively’ distinguish object based on their 
ISM properties. 

On the other hand, the slope, 5, is an important pa¬ 
rameter that connec ts the various physical parameters 
in the ISM (see e.g.. iNiklas fc BecBIlQQTt iDumas et al.l 
I 2 OIII : iSchleicher fc Beckll201.1[ l. Note that at the heart 
of the radio-FIR correlation lies the interdependence 
between various ISM parameters and not their abso¬ 
lute values, viz., c oupling between B and gas density 
(B oc see e.g., iHelou fc BicavlllQ^ INiklas fc Becld 
I 1997 I: iGroyes et ^ I2003I1 . the Kennicutt-Schmidt law 
n ennicu ttl ll998ll and spectrum of the GREs. Therefore, 
the evolution in slope can throw meaningful insights into 
understanding the cause and/or evolution of the radio- 
FIR correlation. In our study, a single slope is enough to 
fit the radio-FIR correlation throughout the luminosity 
and redshift range. We do not find any evidence of its 
evolution within the uncertainty of our measurements. 


5.4. Heterogeneity in the sample 

In spite of the wealth of observational evidence, a clear 
understanding of the reason behind the radio-FIR cor¬ 
relation, across various galaxy types with wide range 
of star-formation activity, magnetic field strengths, na¬ 
ture of ISM turbulence, etc., remains elusive. One of 
the major areas of progress in recent years has been 
in our understanding of how star-formation/supernova 
driven turbulence helps in amplification of magnetic 
field strength in galaxies (iBreitschwerdt et al.l 2009 : 
de Avillez fc Brcitschwerdf] 120051 iMac Low et al.l 2005; 
Gent et al.ll2018;J lblb and how the magnetic field in turn 
couples with the gas density, and hence the star fo rma- 
tion ra te (iGho fc Vishniadl2000l:j^Groves et al.ll200'^ . Re¬ 
cently, ISchleicher fc BeckI (I20I3D interpreted the radio- 
FIR correlation as a result of turbulent magnetic field 
amplification driven by star-formation activity and their 


interplay. Cosmic ray energy loss mechanisms and mag¬ 
netic field strength evolution at higher redshifts may af¬ 
fect the synchrotron emission and hence change the non- 
thermal spectral index giving rise to modification of the 
radio-FIR correlation. 

We would like to point out, since our sample is 
flux-limited in the optical, we have a wide range of 
galaxies in our sample. In terms of Txir, galax¬ 
ies are classified as, 1 ) normal star-forming (Txir ^ 
IQIILq), 2) LIRGs < Txir < 10^^), 3) 

ULIRGs (IO^^Lq < Lxir < 10^^) and 4) Hyper 
LIRGs (Txir > 10^^) (see e.g., iSanders fc Mirabell 
[19^ . As per these definitions, the stacked sources 
mostly fall under the category of normal star-forming 
galaxies. Although, the number of (U)LIRGs are sig¬ 
nificantly less in the local universe (z < 0.3), they 
could contribute sig nificantly to the infrared luminosity 
at hi g her redshifts (iKim fc SandersI 119981: ICaouti et al.l 
120071 : iKilerci Eser et al.1 1201411 . IMagnelli et al.l (|2009f) 
found significant evolution of the contribution from 
(U)LIRGs to the cosmic star formation rate density with 
respect to that from main sequence galaxies. Thus, to¬ 
wards the higher redshift bins (z > 0.9) there could be 
significant contribution from LIRGs and ULIRGs to the 
stacked luminosity. On the other hand, the majority 
of directly detected sources falls under the LIRG and 
ULIRG category. 

In this study, the stacked blue galaxies and the directly 
detected luminous galaxies are expected to have different 
nature of ISM turbulence and magnetic field amplifica¬ 
tion. For normal star-forming galaxies, at lower lumi¬ 
nosity range, the turbulence is mostly driven by star- 
formation activity while for the luminous galaxies the 
tur bulence could be dr i ven by interactions and merg¬ 
ers (jwnii^)iE^[200l; [KHircr^^FiLinilQll)- Even 
perhaps, the mechanisms driving the radio-FIR corre- 
lation can differ in these wide variety of galaxies (see 
ILacki et aI1l2010ll . Ideally, one should study the different 
population of galaxies based on a well defined luminosity 
or stellar-mass selection. Unfortunately, due to the lack 
of sufficient number of directly detected galaxies of sim¬ 
ilar luminosities/stellar-masses across all redshift bins, 
our sample does not allow us to study this aspect. But, 
nonetheless, notable fact is that both class of galaxies fol¬ 
low the radio-FIR correlation with similar parameters. 
This brings to light the universal nature of the radio-FIR 
correlation and suggests similar global processes at play 
in galaxies with very different ISM properties. 

The stacking technique described in this paper trades 
off sample size for improved signal to noise. Such a com¬ 
promise has the obvious drawback that one can only 
study the mean properties of a sample. Outliers in the 
parameter space that correspond to more extreme phys¬ 
ical environments are completely missed out. 

In order to study the radio-FIR correlation in nor¬ 
mal galaxies at high redshift without the above trade¬ 
off requires large, spectroscopically confirmed samples 
of such galaxies spanning the entire redshift range of 
interest. Multi-object spectrographs on large opti¬ 
cal telescopes are now routinely producing such sam- 
ples totalling h undreds of thousands of gal axies (e.g., 
Coil^.e^^ ][2irn] [PRIM US]; IGuzzo et a,l.l[ 2 nil fVTPERS]: 
Newman et al.l 120131 [DEEP2]; iLe Fevre et al.l 120051 







































































Radio-FIR correlation 


13 


[VVDS]; iBaldrv et alJ l2010l [GAMA]). Radio followup 
at sufficient d epth will also be routine with the JVLA 
(| Jarvis et al.l 1201411 and deep surveys with up com- 
ing facilities such as ASKAP l)Norris et al.l 1201 Ih and 
MeerKAT. The situation with deep FIR surveys is less 
hopeful. The Level 5 and Level 6 HerMES data are al¬ 
ready sufficiently confused that increased depth in the 
FIR can only be achieved with a larger telescope oper¬ 
ating at somewhat shorter wavelengths. Such a super- 
Spitzer will need to have a large aperture and instru¬ 
ments far more sensitive than MIPS. Unfortunately, such 
a telescope is nowhere on the horizon. 

The stacking methods described in this paper can 
be gainfully applied to a number of other deep fields 
where almost all the required data are available in public 


archives. 
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APPENDIX 

STACKING ANALYSIS 

Results of our stacking analysis for various wavebands are given in Table O Column 2 shows the range of redshift 
bin and column 3 shows the number of PRIMUS blue galaxies in each bin. In column 4 we list the mean redshift of 
the galaxies. Columns 5 and 6 lists the mean stacked flux density and the rms noise of the stacked images in mJy, 
respectively. In column 7 we list the luminosity derived from stacking in the luminosity space in units of 10^^ W Hz“^. 
For comparison, we list the luminosity derived from stacked flux density and the mean redshift in column 8. We have 
used the values from column 7 in our analysis. The errors on mean stacked flux density and luminosity are the la 
bootstrap errors. 


SED FITS TO STACKED IMAGES 

We show the two components of the SED fitting in the FIR regime. In Figure HI the greybody component is shown 
as a dashed red line and the mid-IR power law component is shown as a blue dashed-dot line. The solid black curve 
shows the total fit. 
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TABLE 3 

Results of stacking analysis. 


Survey 

Redshift 

-^obj 

2^mean 

Mean flux 

Stack image 

Luminosity (xlO^^ W Hz 


bin 



density (mjy) 

rms (mJy) 

Stacked 

Mean 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

0.325 GHz GMRT 

0-0.1 

1525 

0.064 

0.027±0.020 

0.0088 

0.06±0.04 

0.03±0.02 


0.1-0.2 

2378 

0.155 

0.024±0.006 

0.0052 

0.18±0.04 

0.16±0.04 


0.2-0.3 

2836 

0.252 

0.037±0.007 

0.0055 

0.71±0.15 

0.74±0.13 


0.3-0.4 

3268 

0.349 

0.029±0.002 

0.0058 

1.24±0.18 

1.23±0.10 


0.4-0.5 

4410 

0.451 

0.031±0.004 

0.0042 

2.37±0.29 

2.36±0.31 


0.5-0.6 

4476 

0.552 

0.030±0.005 

0.0038 

3.88±0.73 

3.74±0.58 


0.6-0.7 

5339 

0.646 

0.030±0.005 

0.0041 

5.53±0.56 

5.50±0.97 


0.7-0.8 

4374 

0.749 

0.034±0.005 

0.0043 

8.79±0.85 

8.86±1.23 


0.8-0.9 

3548 

0.848 

0.035±0.009 

0.0035 

12.41±3.38 

12.54±3.26 


0.9-1.0 

2114 

0.944 

0.046±0.013 

0.0059 

20.98±4.96 

21.46±6.29 


1.0-1.1 

1348 

1.047 

0.046±0.013 

0.0067 

27.21±9.04 

27.68±8.11 


1.1-1.2 

1160 

1.145 

0.046±0.011 

0.0072 

36.10±10.02 

35.09±8.71 

1.4 GHz VLA 

0-0.1 

722 

0.064 

0.017±0.005 

0.0026 

O.OlitO.OO 

0.02±0.00 


0.1-0.2 

1251 

0.156 

0.012±0.005 

0.0018 

0.08±0.04 

0.08±0.03 


0.2-0.3 

1384 

0.250 

0.013±0.002 

0.0019 

0.27±0.07 

0.25±0.04 


0.3-0.4 

1639 

0.350 

O.OlOiO.OOl 

0.0015 

0.43±0.05 

0.43±0.06 


0.4-0.5 

2450 

0.451 

0.008±0.003 

0.0014 

0.65±0.09 

0.64±0.22 


0.5-0.6 

2664 

0.552 

0.009±0.002 

0.0013 

1.12±0.33 

1.13±0.19 


0.6-0.7 

3040 

0.647 

0.012±0.003 

0.0012 

2.25±0.59 

2.25±0.47 


0.7-0.8 

2270 

0.749 

0.009±0.009 

0.0029 

2.26±1.97 

2.27±2.37 


0.8-0.9 

1941 

0.848 

0.011±0.002 

0.0014 

4.14±0.76 

4.12±0.65 


0.9-1.0 

1079 

0.943 

0.011±0.002 

0.0028 

5.16±1.54 

5.11±1.17 


1.0-1.1 

692 

1.048 

0.012±0.003 

0.0029 

7.67±1.86 

7.54±1.66 


1.1-1.2 

579 

1.145 

0.010±0.006 

0.0029 

7.93±3.93 

7.92±4.32 

24/xm SWIRE 

0-0.1 

1524 

0.064 

0.027±0.004 

0.0016 

0.04±0.005 

0.03±0.004 


0.1-0.2 

2370 

0.155 

0.062±0.004 

0.0012 

0.46±0.03 

0.41±0.02 


0.2-0.3 

2829 

0.252 

0.077±0.004 

0.0013 

1.57±0.10 

1.51±0.08 


0.3-0.4 

3263 

0.349 

0.070±0.002 

0.0013 

2.88±0.16 

2.93±0.08 


0.4-0.5 

4403 

0.451 

0.054±0.001 

0.0012 

4.12±0.16 

4.16±0.11 


0.5-0.6 

4471 

0.552 

0.050±0.002 

0.0010 

6.23±0.36 

6.22±0.25 


0.6-0.7 

5325 

0.646 

0.048±0.002 

0.0009 

8.92±0.38 

8.85±0.34 


0.7-0.8 

4352 

0.749 

0.058±0.002 

0.0012 

15.32±0.51 

15.31±0.43 


0.8-0.9 

3539 

0.848 

0.056±0.002 

0.0010 

20.27±0.56 

20.22±0.68 


0.9-1.0 

2109 

0.944 

0.064±0.001 

0.0021 

30.32±1.04 

30.16±0.69 


1.0-1.1 

1345 

1.047 

0.059±0.003 

0.0023 

35.86±1.73 

35.91±2.04 


1.1-1.2 

1158 

1.145 

0.044±0.003 

0.0022 

33.57±1.93 

33.55±2.48 

TO/im SWIRE 

0-0.1 

1519 

0.064 

0.469±0.080 

0.0963 

0.64±0.12 

0.46±0.08 


0.1-0.2 

2365 

0.155 

0.913±0.058 

0.0871 

6.57±0.49 

6.00±0.38 


0.2-0.3 

2824 

0.252 

1.028±0.088 

0.0709 

19.87±1.76 

20.19±1.72 


0.3-0.4 

3257 

0.349 

0.821±0.053 

0.0553 

34.10±2.15 

34.23±2.22 


0.4-0.5 

4394 

0.451 

0.658±0.056 

0.0580 

49.87±2.79 

50.41±4.25 


0.5-0.6 

4457 

0.552 

0.571±0.034 

0.0600 

71.30±6.20 

71.28±4.27 


0.6-0.7 

5321 

0.646 

0.538±0.067 

0.0517 

98.21±5.50 

98.50±12.17 


0.7-0.8 

4340 

0.749 

0.576±0.040 

0.0672 

152.96±13.28 

152.22±10.51 


0.8-0.9 

3534 

0.848 

0.401±0.058 

0.0538 

143.00±15.29 

144.13±20.83 


0.9-1.0 

2103 

0.944 

0.809±0.308 

0.0784 

377.15±131.50 

379.72±144.33 


1.0-1.1 

1344 

1.047 

0.496±0.094 

0.0765 

307.86±38.37 

301.52±57.37 


1.1-1.2 

1156 

1.145 

0.309±0.046 

0.1277 

232.92±43.61 

234.38±35.20 

160;/m SWIRE 

0-0.1 

1501 

0.064 

1.998±0.387 

0.5270 

2.88±0.46 

1.95±0.38 


0.1-0.2 

2318 

0.155 

3.099±0.368 

0.4853 

22.81±2.23 

20.30±2.41 


0.2-0.3 

2776 

0.252 

2.694±0.242 

0.3390 

53.82±7.23 

52.88±4.74 


0.3-0.4 

3190 

0.349 

2.981±0.243 

0.3311 

125.05±10.65 

124.26±10.13 


0.4-0.5 

4310 

0.451 

2.224±0.259 

0.2910 

167.28±18.95 

170.39±19.83 


0.5-0.6 

4392 

0.552 

2.083±0.199 

0.3003 

261.33±24.83 

260.11±24.90 


0.6-0.7 

5223 

0.646 

2.003±0.194 

0.2349 

374.07±38.55 

366.41±35.47 


0.7-0.8 

4268 

0.749 

1.957±0.198 

0.3171 

526.54±46.87 

517.24±52.21 


0.8-0.9 

3489 

0.848 

2.103±0.195 

0.4453 

761.58±91.49 

754.85±69.84 


0.9-1.0 

2074 

0.944 

2.165±0.567 

0.5286 

1020.48±302.42 

1016.19±266.30 


1.0-1.1 

1316 

1.048 

2.194±0.321 

0.6136 

1372.73±298.74 

1334.89±195.16 


1.1-1.2 

1142 

1.144 

1.946±0.315 

0.5715 

1474.78±371.85 

1477.01±239.37 

250/J,m HerMES 

0-0.1 

1525 

0.064 

0.634±0.139 

0.2008 

1.18±0.22 

0.62±0.14 


0.1-0.2 

2378 

0.155 

2.205±0.183 

0.1671 

17.29±1.43 

14.49±1.21 


0.2-0.3 

2836 

0.252 

2.422±0.224 

0.1579 

48.00±3.96 

47.60±4.41 


0.3-0.4 

3268 

0.349 

2.691±0.160 

0.1440 

111.75±6.68 

112.17±6.68 


0.4-0.5 

4410 

0.451 

2.431±0.147 

0.1208 

186.92±8.69 

186.22±11.23 


0.5-0.6 

4476 

0.552 

2.348±0.146 

0.1269 

292.77±16.40 

293.08±18.20 


0.6-0.7 

5339 

0.646 

2.347±0.086 

0.1219 

433.28±21.38 

429.42±15.68 


0.7-0.8 

4374 

0.749 

2.713±0.131 

0.1324 

714.79±34.14 

717.13±34.72 


0.8-0.9 

3548 

0.848 

2.360±0.146 

0.1560 

843.22±38.97 

847.67±52.30 


0.9-1.0 

2114 

0.944 

2.913±0.247 

0.1885 

1375.33±97.87 

1366.98±115.95 


1.0-1.1 

1348 

1.047 

2.742±0.203 

0.2366 

1680.61±117.42 

1667.16±123.22 


1.1-1.2 

1160 

1.145 

2.461±0.213 

0.2277 

1864.39±166.36 

1869.03±161.61 
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TABLE 3 

( continued ...) Results of stacking analysis. 


Survey 

Redshift 

A^obj 

2mean 

Mean flux 

Stack image 

Luminosity (xlO 

W Hz"^) 


bin 



density (mjy) 

rms (mJy) 

Stacked 

Mean 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

350/xm HerMES 

0-0.1 

1525 

0.064 

0.471±0.126 

0.2139 

0.77±0.21 

0.46±0.12 


0.1-0.2 

2378 

0.155 

1.038±0.155 

0.1972 

8.02±0.89 

6.82±1.02 


0.2-0.3 

2836 

0.252 

0.818±0.101 

0.1886 

15.64±1.90 

16.08±1.99 


0.3-0.4 

3268 

0.349 

1.141±0.135 

0.1475 

45.63±7.68 

47.56±5.62 


0.4-0.5 

4410 

0.451 

1.257±0.098 

0.1158 

96.42±8.35 

96.30±7.50 


0.5-0.6 

4476 

0.552 

1.233±0.080 

0.1266 

155.87±15.47 

153.91±9.96 


0.6-0.7 

5339 

0.646 

1.470±0.053 

0.1269 

272.28±21.03 

269.10±9.71 


0.7-0.8 

4374 

0.749 

1.821±0.145 

0.1423 

484.65±40.00 

481.47±38.26 


0.8-0.9 

3548 

0.848 

1.781±0.156 

0.1591 

635.81±46.56 

639.66±56.10 


0.9-1.0 

2114 

0.944 

1.999±0.224 

0.1679 

948.50±73.44 

938.25±105.05 


1.0-1.1 

1348 

1.047 

2.386±0.160 

0.2496 

1455.06±119.08 

1450.73±97.26 


1.1-1.2 

1160 

1.145 

2.278±0.268 

0.2538 

1730.19±190.89 

1729.98±203.59 

500/J.m HerMES 

0-0.1 

1525 

0.064 

0.599±0.134 

0.2326 

0.87±0.16 

0.59±0.13 


0.1-0.2 

2378 

0.155 

0.368±0.131 

0.1971 

2.88±0.52 

2.42±0.86 


0.2-0.3 

2836 

0.252 

0.380±0.081 

0.1776 

9.11±1.50 

7.46±1.60 


0.3-0.4 

3268 

0.349 

0.275±0.099 

0.1800 

13.02±4.35 

11.46±4.11 


0.4-0.5 

4410 

0.451 

0.451±0.070 

0.1235 

37.97±9.63 

34.53±5.37 


0.5-0.6 

4476 

0.552 

0.474±0.099 

0.1356 

59.29±11.51 

59.15±12.41 


0.6-0.7 

5339 

0.646 

0.809±0.080 

0.1234 

151.81±23.11 

148.13±14.64 


0.7-0.8 

4374 

0.749 

0.941±0.121 

0.1507 

249.93±40.49 

248.61±32.10 


0.8-0.9 

3548 

0.848 

0.850±0.108 

0.1703 

298.76±41.09 

305.34±38.78 


0.9-1.0 

2114 

0.944 

1.084±0.192 

0.1843 

518.44±71.76 

508.75±89.87 


1.0-1.1 

1348 

1.047 

1.288±0.159 

0.2314 

789.56±122.30 

783.08±96.78 


1.1-1.2 

1160 

1.145 

1.268±0.151 

0.2526 

960.91±178.99 

963.30±114.79 
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Fig. 8.— Far-infrared SED fit to the stacked images in various redshift bins. The red dashed curve represents the modified black body 
component fit and the blue dash-dot line represents the mid-infrared power law fit. The black solid line represents the sum of these 
components and provides an excellent fit to the data in almost every case. 






























